Neutrinos from Gamma-Ray Bursts in Pulsar Wind Bubbles: ~ 10 16 eV 
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The supranova model for 7-ray bursts (GRBs) is becoming increasingly more popular. In this 
scenario the GRB occurs weeks to years after a supernova explosion, and is located inside a pulsar 
wind bubble (PWB). Protons accelerated in the internal shocks that emit the GRB may interact 
with the external PWB photons producing pions which decay into ~ 10 16 eV neutrinos. A km 2 
neutrino detector would observe several events per year correlated with the GRBs. 
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The leading models for Gamma-Ray Bursts (GRBs) 
involve a relativistic outflow emanating from a compact 
central source. The ultimate energy source is rapid accre- 
tion onto a newly formed stellar mass black hole. Long 
duration (> 2 s) GRBs, which include all GRBs with 
observed afterglow emission and ~ 2/3 of all GRBs, are 
widely assumed to originate from a massive star progen- 
itor. This is supported by growing evidence for GRBs 
occurring in star forming regions within their host galax- 
ies Q . The leading model for long duration GRBs is the 
collapsar model |2(, where a massive star promptly col- 
lapses into a black hole, and forms a relativistic jet that 
penetrates through the stellar envelope and produces the 
GRB. An interesting alternative (though somewhat more 
debated) model is the supranova model 0, where a su- 
pernova explosion leaves behind a supra-massive neutron 
star (SMNS), of mass ~ 2.5 — 3 M®, which loses its rota- 
tional energy on a time scale, t s( j ~ weeks to years, and 
collapses to a black hole, triggering the GRB. The most 
natural mechanism by which the SMNS can lose its rota- 
tional energy is through a strong pulsar type wind, which 
is expected to create a pulsar wind bubble (PWB) 0,O|. 

The prompt 7-ray emission in GRBs is believed to arise 
from internal shocks within the relativistic outflow, of av- 
erage isotropic luminosity L ~ 10 52 erg/s, that form due 
to variability in its Lorentz factor, T, on a time scale 
t v . These shocks accelerate electrons to relativistic ran- 
dom velocities, which then emit synchrotron (and per- 
haps also synchrotron sclf-Compton) radiation, that con- 
stitutes the observed 7-ray emission. In the region where 
electrons are accelerated, protons are also expected to be 
shock accelerated: the conditions in the dissipation re- 
gion allow proton acceleration up to £ p , ma x ~ 10 20 eV 0. 

In this Letter we consider the neutrino production in 
GRBs that occur inside PWBs, as is expected in the 
supranova model for GRBs. The neutrino production via 
photomeson interactions between the relativistic protons 
accelerated in the internal shocks and the synchrotron 
photons that are emitted in these shocks, has already 
been calculated || , while p-p collisions were shown to be 
relatively unimportant. Therefore, we focus on neutrino 
production through photomeson interactions with the ex- 
ternal photons from the rich radiation field inside the 
PWB. This process turns out to be dominant for a wide 



range of parameters: t s( j < 0.2 yr for a typical GRB, and 
i s d < 2 yr for X-ray flashes. The neutrinos produced via 
this mechanism have energies e v ~ 10 15 — 10 17 (10 19 ) eV 
for typical GRBs (X-ray flashes) and are emitted simul- 
taneously with the prompt 7-ray (X-ray) emission. Their 
energy spectrum consists of several power law segments, 
and its overall shape depends on the model parameters. 
For t s d < 0.1 yr, the v's would not be accompanied by 
a detectable GRB. We find that the ~ 10 16 eV neutrino 
fluence from a GRB at z ~ 0.1 — 1 implies ~ 0.01 — 1 up- 
ward moving muons in a km 2 detector simultaneous with 
the 7-rays. The neutrino signal from an individual GRB 
that is pointed towards us, and is therefore detectable in 
7-rays, is above the atmospheric neutrino background. 

Recently, calculated the neutrino emission in the 
supranova model from protons that escape the internal 
shocks region and produce pions inside the SNR shell. 

Neutrino production in the GRB. - The internal 
shocks in GRBs are believed to be mildly relativistic. 
Therefore, the proton energy distribution should be close 
to that for Fermi acceleration in a Newtonian shock, 
dnp/dsp oc e~ 2 . Moreover, the power law index of the 
electron and proton energy spectra should be the same, 
and the values inferred for the electrons (from the prompt 
GRB spectrum) are dn e /de e oc e~ p with p ~ 2 — 2.5. 



We therefore adopt dn p /de p 
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quantities are measured in the comoving (lab) frame. 
Protons of energy e p interact mostly with photons that 
satisfy the A resonance condition, e p .a — 0.3GeV 2 /e 7 . 
The minimal photon energy relevant for P7 interactions 
is £ 7 ,min = 0.3 GeV 2 /e Pimax ~3x 10~ 3 eV. For reason- 
able parameters [i s d,-i = £sd/(0.1yr) > 0.26] this falls 
above the self absorption frequency |5J , and the relevant 
part of the PWB spectrum consists of two power laws: 
dn-f/de-f oc e 7 3 ^ 2 for e 1 < e 7 & = hvb m and oc e 7 s ^ 2 1 for 

e 7 > e 7 b, where Uf, m ~ 1.6 X lO 15 ^ 3 ^ Hz is the peak 
frequency of the PWB vF v spectrum, n 7 is the photon 
number density, and s ~ 2.2 is the power law index of 
the PWB electrons. Photons of energy e 7 b satisfy the A 
resonance condition with protons of energy 
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where (3 b c = O.l/^-ic is the velocity of the SNR shell, 
-f w = 10 4 - 5 7 W! 4.5 and E mt — 10 53 -E 53 erg are the Lorentz 
factor and total energy of the pulsar wind, £ e = £ e /3/3 
(r) = 772/32/3) is the fraction of the wind energy in 
pairs (protons), and e be = e 6e / 3 /3 (e bB = 10~ 3 e 6B _ 3 ) is 
the fraction of P WB energy in electrons (magnetic field) . 
The corresponding neutrino energy is e vb « e p h/20 ~ 
2 x lO 15 ^ 2 —1 eV, similar to that expected from the P7 
interactions with the GRB photons 0, for tj sc i ~ 0.1 yr. 
The GRB emission can be detected simultaneously with 
these neutrinos only if the Thompson optical depth is 
< 1, which can be obtained for t & & ~ 0.1 yr only if the 
SNR shell is clumpy 5j. In the case of a uniform shell 
we need 7j s( j > 0.4 yr, and in turn e vb > 2 x 10 16 eV. 

The internal shocks occur over a distance AR ~ R = 
2T 2 ct v . Thus, the optical depth to photo-pion production 
at the A resonance, for protons of energy e p , is 
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where r 2 . 5 = T/10 2 5 , 

e 7 = 0.3 Gev 2 /e p , / = / 1/3 /3 is the fraction of the PWB 
radius up to which most of its radiation is emitted, and 
(3 = s/2 (1/2) for e p < e pb (e p > e pb ) is the spectral slope 
of the seed PWB synchrotron photons. The fraction of 
the proton energy that is lost to pion production is |10| 



jW(£p) ~ 1-exp [-r P7 (ep)/5] 



,(e p )/5] . (3) 



The factor 5 is since the proton loses ~ 0.2 of its energy in 
a single interaction. We denote e p for which f VK {e p ) « 1 
by s pT = 10 18 e pT i8 eV [i.e. T pl {e pT ) = 5], and obtain 
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decay of charged pions created in interactions be- 
PWB photons and GRB protons, produces high 
neutrinos, 7r + — * ji + + — > e + + v e + Vfi + v^, 
each neutrino receives ~ 5% of the proton energy, 
energy of the protons accelerated in the internal 
of GRBs is expected to be similar to the 7-ray 
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FIG. 1: The proton energies that correspond to break ener- 
gies in the neutrino spectrum, e„ ~ e p /20, as a function of 

tad- 



decays. We have taken into account that although the 
initial neutrino flavor ratio $ Ue : : is 1 : 2 : 0, 
neutrino oscillations cause it to be 1 : 1 : 1 at the Earth. 

The pions may lose some energy via synchrotron or 
inverse- Compton (IC) emission. If these energy losses 
are important, then the final energy left in the ir + when 
it decays, out of which 3/4 goes to neutrinos, can be much 
smaller than its original energy. We find that f 7TV {e p ) rs 
1 - exp(-t; ad /r;) » min(l,4 d /T;), where < w 2.6 x 
10 8 e' 7T /(m 7T c 2 ) s is the life time of the pion, and t' Tad = 
(tgyn + tjQ ) _1 ~ mm(t' syn , t[ c ) is the time for radiative 
losses due to both synchrotron and IC losses, the time 
for which are t' syn and t' IC , respectively. We have 
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where U' B = B 



7r is the energy density in the mag- 
netic field, within the shocked fluid behind the internal 
shocks, and U'^e'^) is the energy density of photons be- 
low the Klein-Nishina limit, e' K N — (m.c 2 ) 2 /^ IC 
losses due to scattering of the GRB photons were shown 
to be unimportant We consider the IC losses due to 
the upscattering of the external PWB photons, and find 
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where e e and e# are the equipartition parameters of the 
GRB, and e^ig = e^/lO 18 eV. The radiative losses be- 
come important for t' r&d < t' v , which corresponds to 

P ^ £ p s 
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energy in 7-rays, f w {e p ) = f pv (e p )f^{e p ), while f p ^(e p ) e «yn g _ 2 3 e i/2 £ -i/2 i -i/2 r 4 ^ ^ 
is given in Eq. and f.Ki,{e p ) is the fraction of the 
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FIG. 2: The muon neutrino spectrum, for our fiducial param- 
eters and t Bd = 0.01, 0.07, 0.4, 30 yr, for (F,t v ,- 2 ) = (10 2 5 , 1) 
(solid), which correspond to the 4 different orderings of the 
break energies, and (T,t V] -2) = (600,5) (dashed). For com- 
parison we show the spectrum of for (r, ^,-2) = (10 2,5 , 1) 
(dot-dashed) and (600, 5) (dotted), where the thick (thin) lines 
are with (without) the effects of the PWB radiation field (that 
inflict energy losses on the protons, pions and muons). 



The protons may also lose energy via P7 interactions with 
the GRB photons [8j, However, r p7 for this process is 
typically < 1, so that it does not have a large effect on P7 
interactions with the PWB photons, on which we focus. 

Since the life time of the muons is ~ 100 times longer 
than that of the pions, they would have significant radia- 
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FIG. 3: The number of muon events per GRB in a km 2 
detector, as a function of t s <x, for or fiducial parameters and 
three different values of (T,t v ). The horizontal lines show 
the events expected due to P7 interactions with the GRB 
photons Q. 



fixed value of t s d,-i = lj ^2.5 ^ 1 implies r p7 > 1 and 
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const ~ 10 eV, while T 2 .5 < 1 implies 
E ps = £p^ n oc r 4 . This implies a stronger neutrino emis- 
sion which can reach higher energies for larger values of 
r. Since a larger T implies a lower typical synchrotron 
frequency for the prompt GRB, this may apply to X- 
ray flashes, if they are indeed GRBs with relatively large 
Lorentz factors T and/or a large t v [TT| . 

As can be seen from Fig. ^ for relevant parameters, 
there are four different orderings of these break energies: 

(i) £„ T < £ v b < E vs , (ii) E v b < E VT < E vs , (Hi) E v b < E vs < 

E vr , (iv) e vs < E v b < e vt . Each of these orderings results 
in a different shape for the spectrum, that consists of 3 
or 4 power laws, as can be seen in Fig. 2 (solid line). The 
peak of the £ 2 (dN v / d,E v ) spectrum, in units of /q, is 
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tive losses at an energy of Ens ~ £*- s /10 « £ ps /50. This s(A-s) , , A s / 2 Jr u/r f <r r t <r f 



causes a reduction of up to a factor of 3 in the total neu- 
trino flux in the range ~ (0.1 — l)£i/ S , since only that 
are produced directly in 7r + decay contribute consider- 
ably to the neutrino flux. Note that since both ratios in 
Eqs. (0 and © scale as £~ 2 , we always have t' lad /r^ oc 
£~ 2 , and therefore the spectrum steepens by a factor of 
(e v /e us )~ 2 for e v > e vs . This is also evident from the 
fact that f„ v (e p « 5^) « min(l, t' syn /T^, t[ c /r^). 

Fig. ^ shows the proton energies corresponding to the 
neutrino break energies e vs , e v \, and £„ r s=s £ pr /20, as 
a function of t s &. From Eqs. ©, © and (JTDJl, we see 

that 



For example, spectrum (iii) applies for i s d = 0.4 yr, where 
the GRB is detectable, and is given by 
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Fig. [5]shows the muon neutrino spectrum for our fiducial 
parameters (solid line). The spectrum of the other neu- 
trino flavors is the same, while the spectrum of the 7-rays 
from the pion decay is almost the same, just substitut- 
ing e 7 w 2e u and with a normalization larger by a factor 
of 4. In Fig. 12 we also compare our neutrino spectra 
with that of y| for different internal shock parameters. 
We see that for reasonably large values of T and t v , the 
PWB photons are the dominant target for photomeson 
interactions for £ s( j < 4 yr. The afterglow neutrino emis- 
sion with e„ > 10 17 eV 0] should be distinguishable 
from the prompt emission (from P7 interaction with ei- 
ther GRB or PWB photons) due to the very different 
spectrum, and may also be delayed in time. 

Implications — The high energy, ~ 10 16 eV, neutrinos 
produced by photomeson interactions between protons 
accelerated in the internal shocks of GRBs and the PWB 
photons may be detected by the future neutrino tele- 
scopes 0]. The probability, P Ufl , that a neutrino would 
produce a high energy muon in the detector is approxi- 
mately given by the ratio of the high energy muon range 
to the neutrino mean free path. For the neutrinos stud- 
ied here, P„ M w 1.3 x lfT^/lO 3 TeV)' 3 , withj = 1 for 
e v < 10 3 TeV and 0=1/2 for e v > 10 3 TeV 

In Fig. we report the expected number of events 
in a km 2 detector as a function of t s( j, for our fiducial 
parameters, for a GRB at z ~ 1 with i? 7 ,i SO = 1 53 ergs. 
We consider three different sets of values for V and t v , 
and compare the resulting numbers with the expected 
number of events due to photomeson interactions with 
the GRB photons 0. For a typical GRB with T = 10 25 
and t v — 10 ms, v's from interactions with PWB photons 
dominate over v's from GRB photons for i S( j < 0.2 yr. 

GRBs that occur inside spherical PWBs with t s( j ~ 
0.1 — 0.2 yr would have a peculiar and short lived after- 
glow emission Since they occur for a wide range of 
i s d values, we expect their rate to be similar to that of 
typical GRBs (i.e. ~ 10 3 yr -1 beamed toward us). For 
these values of t s d we expect ~ 0.01 — 0.1 events per burst 
corresponding to the detection of several tens of neutrino 
induced muons per year. These neutrino bursts should 



be easily detected above the background, since the neu- 
trinos would be correlated, both in time and direction, 
with the GRB 7-rays. Note, that at the high energy con- 
sidered, knowledge of burst direction and time will allow 
to discriminate the neutrino signal from the background 
by looking not only for upward moving neutrino induced 
muons, but also by looking for down-going muons. For 
larger values of 0.2 yr < t s d < 1 yr, the neutrino flux 
due to photomeson interaction with the PWB photons 
will dominate over the one due photomeson interaction 
with the GRB photons @, if the GRB has T > 600 and 
tv ^ 50 ms. Since larger T and t v imply a lower typi- 
cal synchrotron frequency for the prompt GRB, this may 
apply to X-ray flashes, if they are indeed GRBs with rel- 
atively large Lorentz factors and/or a large variability 
time, t v For these events no afterglow emission has 
been detected and this can be explained considering the 
fact that for i s d < 1 yr the GRB would have a peculiar 
and short lived afterglow emission. The typical neutrino 
energy is expected to be in the range ~ 10 15 — 10 17 eV 
for typical GRBs, and - 10 15 - 10 19 eV for X-ray flashes. 

We expect 10~ 5 — 0.01 events per X-ray flash for this 
range of i s d corresponding to a detection of 0.01 — 10 
events per year. Again, these neutrino bursts should be 
easily detected above the background, since the neutrinos 
would be correlated, both in time and direction, with the 
X-rays of the X-ray flashes. This neutrino emission would 
be simultaneous with the 7-ray/X-ray emission from the 
GRB and should be easily distinguishable from neutrinos 
emitted after the 7-ray phase of the GRB 13] or < 100 s 
before the GRB 141. Detection of high energy neutrinos 
will test the shock acceleration mechanism and the sug- 
gestion that GRBs are the sources of ultra-high energy 
protons, since > 10 16 eV neutrino production requires 
protons of energy > 10 18 eV, and will help to establish 
whether indeed most GRBs occur inside PWBs. 
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